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Single crystal neutron diffraction is combined with synchrotron x-ray scattering to identify the 
different super lattice phases present in Cso. 8 Eei. 6 AS'e 2 . A combination of single crystal refinements 
and first principles modelling are used to provide structural solutions for the \/5 x and a/2 x 
super lattice phases. The \/5 x a/ 5 superlattice structure is predominantly composed of ordered 
Fe vacancies and Fe distortions, whereas the a/2 x a/ 2 super lattice is composed of ordered Cs 
vacancies. The Cs vacancies only order within the plane, causing Bragg rods in reciprocal space. 

By mapping x-ray diffraction measurements with narrow spatial resolution over the surface of the 
sample, the structural domain pattern was determined, consistent with the notion of a majority 
antiferromagnetic \/5 x \/5 phase and a superconducting \/2 x \/2 phase. 


I. INTRODUCTION 

The alkaline iron selenides AxFe 2 -ySe 2 (A = K, Rb, 
Gs) have attracted much interest recently due to the 
observation of superconductivity with Tc in 

conjunction with antiferromagnetism with an unusually 
high ordering temperature Tat of up to 559K and large 
ordered moment of about 3.3/iB per F^. This observa¬ 
tion suggests the possibility of coexistence between these 
two orthogonal phenomena, superconductivity and mag¬ 
netic ordering. Other iron pnictide superconductors have 
phase diagrams that indicate coexistence of magnetism 
and superconductivity in certain regions, however the 
highest Tc’s and bulk superconductivity are only found 
when the magnetic state has been suppressecP®. The 
basic principles of superconductivity indicate that mag¬ 
netic fields cannot permeate a superconducting region. 

At elevated temperatures these compounds exhibit the 
tetragonal space group lA/mmm. Fe vacancies order 
below ~ 600K into a a/5 x ^/b structure with lAjm 
symmetr}®. The Fe lattice in this structure magnetically 
orders into a block antiferromagnet, with neighbouring 
groups of four Fe ions having opposite spin^. Other 
phases have been observed in the system by diffraction 
and scanning electron micros copy, indicating the pres¬ 
ence of a \/2 X ^/2 phasd^^^ and even a 2 x 4 rhombus 
phas^^. This has led to the idea that the superconduc¬ 
tivity in this system may be in a different phase to the 
principal magnetic phase. 

The observation of multiple phases has driven a num¬ 
ber of reports on phase-separation in this system, us¬ 
ing various techniques to probe either the surface or 
the bulk material. Transmission electron microscopy 
(TEM) measurements observed parallel lamellae of or¬ 
dered and disordered Ee-vacancy regionJ^ and syn¬ 


chrotron x-ray diffraction measurements determined the 
nanoscale phase-separation of a magnetic, Ee-vacancy su¬ 
perlattice and a nonmagnetic phase with an in-plane 
compressed lattic^^^l^ Scanning Tunneling Microscope 
(STM) measurements showed that the superconducting 
state would not survive in regions with Ee vacanciesl^ 
and finally bulk measurements with NRM concluded 
that the superconducting and antiferromagnetic proper¬ 
ties are part of two different phase 

In this paper we will describe the results of single crys¬ 
tal diffraction measurements on Cso,^Fei,QSe 2 leading 
to the determination of the a/2 x a/2 superstructure and 
show that this phase is spatially separated from the mag¬ 
netic \/5 X \/5 structure by scanning the sample surface 
with a focused x-ray beam. 


II. METHODS 

Single crystals of CsxFe 2 -ySe 2 were grown by the 
Bridgman process as described in Ref.^^. The nomi¬ 
nal composition of the crystals used in this study is 
Cso_sFei_QSe 2 ^ and their superconducti ng an d magnetic 
properties have been reported previousl}EE2Hini^ sum¬ 
mary, the inclusion of Gs as the intercalation ion leads 
to a large increase in the c-lattice parameter compared 
with KxFe 2 -ySe 2 , despite this however, the reported su¬ 
perconducting transition is roughly the same at Tc ~ 
30K and the antiferromagnetic ordering transition is still 
above room temperature at T/v ~ 500K. The sam¬ 
ples were coated in Gytop varnish to protect them from 
reaction with the atmosphere. A single grain sample 
was measured using a white beam of neutrons on the 
SXD instrument installed at the ISIS spallation neutron 
source (Oxfordshire, UK) for single crystal diffraction 






2 


analysi^. To determine the spatial separation of the 
domains, a second sample from the same crystal growth 
was cleaved under a He atmosphere and measured with 
a focused beam of synchrotron x-rays on beamline 116 at 
Diamond light source (Oxfordshire, UK]P. 

Density Functional Theory (DFT) calculations were 
performed using the plane-wave pseudopotential ap¬ 
proach as implemented in the CASTEP cod^^. The ini¬ 
tial calculations used the block antiferromagnetic struc¬ 
ture previously shown to be the lowest energy state for 
systems with no Fe vacancie^^. A system with two layer 
periodicity was geometry optimised such that the resid¬ 
ual forces were |F| <= 0.036eE/A. This optimised struc¬ 
ture was then used as the base for all further calculations. 
Changes to the Cs layer position were not found to in¬ 
crease the maximum residual force in any layer and con¬ 
sequently no further optimisations were performed. The 
spin polarised Perdew-Burke-Ernzerhof generalised gra¬ 
dient approximation was used to model the exchange and 
correlatiorP^. The default CASTEP ultrasoft pseudopo¬ 
tentials (Version 7) were used with a plane wave cutoff 
of 550 eV. The electronic Brillouin zone was sampled by 
a Monkhorst-Pack gricP^ with a density of at least 0.02 
A. Calculations were performed with a Gaussian smear¬ 
ing of 0.05 eV, giving a total error in the calculation of 
energy of about 0.2 meV. 


III. SINGLE CRYSTAL NEUTRON 
DIFFRACTION 

A. Parent Structures 

A single crystal of CsQ_sFei_QSe 2 was mounted in a low 
background sample environment on SXD and measured 
at the base temperature, 30K. Several orientations were 
measured to build complete coverage of reciprocal space 
and the orientation was chosen so that the principal axes 
were centred in the rear detectors with highest resolution 
in wavevector transfer. The resulting diffraction pattern, 
illustrated in Fig. shows a large number of reflections. 
The larger peaks are from two overlapping tetragonal av¬ 
erage structures and the smaller peaks are satellites pro¬ 
duced by commensurate superstructures. Concentrating 
initially on the principal reflections of the average struc¬ 
tures, the two overlapping phases have slightly different 
lattice parameters. The lattice parameters were deter¬ 
mined by fitting the peak profiles along the principal di¬ 
rections when these directions were aligned at the centre 
of a back scattering detector. Phase 1 has the great¬ 
est intensity and lattice parameters ai = 3.9472(2)A 
and Cl = 15.2401(1) A. The weaker second phase is ex¬ 
panded along the c direction and contracted in a, having 
a2 = 3.8554(1)A and C2 = 15.4931(5)A. The variation of 
the peak intensity indicates the sample fraction of phase 
2 to be 42(3)%. 

These large Bragg peaks were integrated using 3D 
profile fitting methods and refined using the program 




FIG. 1. (Colour online) Single crystal diffraction pattern of 
Cso. 8 Fe 1 . 6 Se 2 measured on SXD. A cut through the (/i,/c,0) 
plane is illustrated on the left in (a) with superlattice peaks 
visible from both phases present in the sample: the \/5 x \/5 
phase illustrated in (b) and the \/2 x y/2 phase in (c). Light 
blue circles show the ^/5 x ^/5 lattice, light red squares the 
\/2 X y/2 lattice. A simulation of the diffraction data on the 
right hand side, using the models in (b) and (c), exhibiting 
good agreement with the measured data. 

Jana200d^. The two phases were indexed, integrated 
and refined independently after absorption corrections 
were applied. Refining the average structure is unable 
to give information on the ordering of vacancies but can 
be used to determine the average occupancies and ther¬ 
mal parameters of different sites. Initially the refine¬ 
ments from both phases gave similar results and were 
consistent with previous literature on this phas^, with 
Cs and Fe concentrations of approximately 80% of the 
fully occupied values. It was realised, however, that a 
number of reflections overlap between the two phases 
which would have an averaging effect on the occupancies. 
At low wavevector transfers, where reflections are closer 
together in reciprocal space, overlapping reflections are 
more likely due to the low detector resolution on SXD. All 
peaks with low wavevector transfer were omitted and a 
cut off of IQ I < 4A“^ was used. The limit of |Q| < 4A“^ 
was chosen as it was possible to distinguish separate pro¬ 
files for overlapping reflections above this limit, as can be 
seen in Fig. Refinements from the resulting reflections 
are presented in Table [Ij 

When low resolution reflections are omitted, the oc- 
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Wavevector Q || a*+b* (A'^) 


FIG. 2. (Colour online) Bragg rods observed in the {h,h,l) 
plane of the neutron diffraction data. The left-hand side gives 
the neutron diffraction and the right-hand side compares the 
simulated data from a model composed of 100 randomly vary¬ 
ing iterations of the Cs-vacancy structure, approximating dis¬ 
order along c. 


cupancies of the two phases become different. The first 
phase refines to a concentration of Cso.97(5)Fei.6o(4)5'e2, 
consistent with a filled Cs site and vacancies on the 
Fe sites. The second phase has a concentration of 
C'so. 63 ( 9 ) 92 ( 8 ) 5'e2, which is consistent with Cs vacan¬ 
cies and filled Fe sites. The increased Fe concentration 
in the second phase is consistent with the increase in 
the c-axis lattice parameter of this phase. R-factors of 
the two solutions were quite high compared with typi¬ 
cal published values at = 13.07% {R = 9.64%) and 
Rw = 24.47% {R = 16.78%) respectively, though their 
large values are attributable to the difficulty in determin¬ 
ing accurate intensities from partially-overlapping reflec¬ 
tions. 

An alternative method of selecting the peaks is to com¬ 
pare the difference in d-spacing between the two compo¬ 
nents of each peak to the resolution function of SXD, 
specified in Fig. 6 of Keen et al^. Omitting reflections 
where the difference in d-spacing between the two sets 
of lattice parameters are less than the resolution at the 
particular 20 value resulted in a similar final set of reflec¬ 
tions to the omission of reflections with |Q| < 4A“^ and 
the refinement results were found to be consistent. 



FIG. 3. (Golour online) Separation of Bragg reflections along 
principal directions in single crystal neutron diffraction data. 
In each case the peak is split into two components; a larger 
component with larger a and smaller c and a weaker compo¬ 
nent. Positions of each set of lattice parameters for phase 1 
and 2 are shown by blue and red dashed lines respectively. In¬ 
set within each plot is an example of a reflection in the same 
direction with |Q| < 4A“^, where the peaks are strongly over¬ 
lapping. 

B. Superstructures 

Weaker peaks were observed away from the principal 
diffraction spots which can be described as commensu¬ 
rate superlattice peaks. The first set of peaks can be 
indexed on a commensurate grid defined by the lattice 
parameters from phase 1 and has propagation vectors: 

^\/5 ~ ~ 

^\/5 = 

^x/5 = 

This phase is named the x ^/b phase, on account 
of the in-plane lattice vectors having a length of ai\/b. 
These propagation vectors determine one of two sym¬ 
metrically equivalent domains, where the second domain 
can be generated by a reflection in either the {1,0,0} 
or {1,1,0} planes. The second set of superlattice peaks 
can only be indexed by lattice parameters of the second 
phase and has propagation vectors: 

^ x /2 = ^2 — ^2 
^\/2 = ^2 + ^2 
%/2 = ^2 

This phase is named the x ^/2 phase, on account of 
the in-plane lattice vectors having a length of a2\/2- The 
superlattice peaks observed in this lattice are not peaked 
out of the plane but were instead found to exhibit rods 
along /, as shown in Fig. This scattering was very 
weak and it was not possible to determine any intensity 
variation along the rod. Such ID scattering occurs from 
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(a) All reflections 



Atom Site Occupancy 

Uii/22 

U33 

Phase 1 

Cs 

2 b 

0.82(6) 

0.036(3) 

0.007(3) 

Refs: 681 

Fe 

Ad 

0.81(1) 

0.0100(6) 0.0033(8) 

= 13.89% 

Se 

Ae 

1 

0.0116(6) 0.0093(8) 

Phase 2 

Cs 

2 b 

0.83(6) 

0.039(6) 

0.021(7) 

Refs: 480 

Fe 

Ad 

0.80(2) 

0 .010(1) 

0.007(2) 

= 23.89% 

Se 

Ae 

1 

0.013(1) 

0.017(2) 

(b) Only reflections with \Q\ > AA ^ 



Atom Site Occupancy 

Ull/22 

U33 

Phase 1 

Cs 

2 b 

0.97(5) 

0.042(4) 

0.014(3) 

Refs: 552 

Fe 

Ad 

0.81(1) 

0 .010(1) 

0.004(1) 

= 13.07% 

Se 

Ae 

1 

0.013(1) 

0 .010(1) 

Phase 2 

Cs 

2 b 

0.63(9) 

0.023(7) 

0 .012(8) 

Refs: 381 

Fe 

Ad 

0.96(4) 

0 .012(2) 

0 .010(2) 

= 24.47% 

Se 

Ae 

1 

0 .011(2) 

0.013(2) 


TABLE L Occupancies and anisotropic thermal parameters 
of the two phases of CsQ,sFei,QSe 2 by refinements of the 
principal Bragg reflections. Both phases are refined in the 
tetragonal space group M/mmm. Phase 1 has lattice pa¬ 
rameters ai = 3.9472(2)A, ci = 15.2401(1)A and phase 2 
02 = 3.8554(1) A and C2 = 15.4931 (5) A. The atomic sites are 
2b = (|, |,0), Ad = (|,0, I), and 4e = (0,0,0.1540(3)). Re¬ 
flections were collected up to a maximum \Q\ = 16A“^, the 
number of reflections observed in each dataset is given by 
“Refs:”. Table (a) gives the refined values when all measured 
reflections are used and table (b) shows the improved values 
given when refined using high resolution reflections. 


a lack of correlation in the direction of the rod, implying 
that successive planes are not ordered with respect to one 
another. 


vice versa. The size of these regions is roughly 300/im x 
300 fim for the dominant first phase and bOfim x 300/im 
for the minority phase. 
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FIG. 4. (Colour online) Scans across the surface of the sam¬ 
ple at the (a) (2, 2,0)i, (b) (2, 2 , 0 ) 2 , (c) (2.6,2.2,0)i and (d) 
(2.5, 2.5,— 1)2 reflections, where the reflections are defined us¬ 
ing lattice parameters from phase 1 or 2. Each element is 
normalised against the sample surface by dividing by the sum 
of the two (2,2,0) reflections. 


IV. SYNCHROTRON X-RAY DIFFRACTION 

Scans over a freshly cleaved surface were performed on 
116 to determine the pattern of structural domains. The 
sample was mounted in reflection geometry in a cryostat 
and measurements were performed at the base tempera¬ 
ture of 8K. The beam size was closed down to 20xl0/im 
and an energy of 7.05keV was chosen to make the x-ray 
penetration depth comparable to the beam size, with an 
attenuation length of ~ 20/im. Scans across the sample 
were performed by aligning on a specific peak and then 
translating the sample through the beam. The Pilatus 
area detector was used to integrate over the peak in 20 
and X, and a rocking scan in 0 was performed at each 
point. Scans were taken at the two components of the 
principal diffraction spot and also at representative su¬ 
perlattice peaks for each phase. The results, shown in 
Fig. i indicate that the regions of greatest intensity in 
phase 1 coincide with the regions of lowest intensity in 
phase 2. The same conclusion is drawn from the super¬ 
lattice peaks, with the ^/b x \/5 peak more intense in 
the regions where the \/2 x \/2 peak is less intense, and 


V. MODELS 

Two separate phases can be distinguished in the exper¬ 
imental data, with different commensurate super lattices. 
The first phase has been well studied previously in the 
literature and can be described by a x ^/5 iron va¬ 
cancy superstructure, illustrated in Figj^b). Simulations 
of this model are consistent with our diffraction data, 
shown in Fig. [^a). The neutron refinements, which sep¬ 
arate the two phases, indicate that the Cs occupancy is 
filled in this phase and indeed, simulating the diffracted 
intensities in each case indicates the relative intensities of 
the superlattice peaks do not change significantly when 
changing the Cs concentration. 

Phase 2 has been observed prev iously i n the litera¬ 
ture as an unknown \/2 x \/2 phasd^^^^^i^. Our refine¬ 
ments indicate that only the Cs site contains vacancies in 
this phase. The commensurate supercell comprises two 
average unit cells, or two occupiable Cs sites per layer. 
The Cs ions sit at high symmetry positions in the parent 
space group, preventing their movement. The only solu¬ 
tion that gives peak positions consistent with the diffrac¬ 
tion pattern is to remove every other Cs ion from the 
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layer as illustrated in Fig. This produces stripes of Cs 
ions and a concentration of 50%, which is consistent with 
the refinements. On consecutive Cs layers, there are two 
possible positions to place the vacancy arising from two 
layout options per layer, see Fig. ^b). It is not possi¬ 
ble to determine, by symmetry, which layout should be 
preferred. Simulations of this model, expanded in the c 
direction over a large number of cells with random ori¬ 
entations of Cs layers produce a diffraction pattern with 
ID diffuse stripes along /, consistent with the neutron 
diffraction data, see Fig. 


(a) (b) 



FIG. 5. (Colour online) Cs layer ordering in the a/2 x a/2 
phase. Image (a) illustrates a single Cs layer, where Cs ions 
are light blue and vacancies are white. They are shown above 
octahedra formed by Fe and Se ions, coloured brown and 
green respectively. The diagram in (b) is a representation 
of the different possible layouts of Cs ions in different layers, 
green and red circles indicate the two possible layouts in the 
first layer and blue and orange squares indicate positions in 
the second layer. The small energy difference between these 
different choices leads to a lack of correlations between each 
layer, producing Bragg rods in the diffraction pattern. 


VI. AB INITIO CALCULATIONS 

Ah initio calculations were performed of the a/2 x a/2 
structure using the CASTEP code. To show that the pro¬ 
posed Cs vacancy structure is the origin of the 2D order 
it is necessary to show two things. First it must be shown 
that the system will order at the very least in-plane and 
second that there is no preferred stacking sequence. Two 
calculations were performed with different in-plane mod¬ 
els, the first with the proposed ordered vacancy struc¬ 
ture shown in Fig. ia) and the second with a single Cs 
displaced to an unfavourable site. To investigate if the 
system will order 3 dimensionally it is then necessary to 
consider stacking sequences. The vacancy placement in 
the first two Cs layers is irrelevant as all four options, AB, 
AD, CD or CB are invariant under either translation or 
inversion symmetry. The third and fourth layers however 
can either sit directly on top of their respective original 
layers or in the previously vacant sites, having a stack¬ 
ing sequence such as ABAB or ABCD. These different 
options are illustrated in Fig. [^b). 


Within an expanded supercell of the Cs superstructure, 
with a' = 4a2, b' = 4^2 and c' = C 2 , a single Cs ion was 
moved to a vacancy site. This movement led to an in¬ 
crease in total energy of 0.851 eV. This is a large increase, 
though most of it is attributable to simple Coulombic re¬ 
pulsion. In the proposed structure each Cs has four near¬ 
est neighbours and a Cs-Cs bond length of 5.36 A. The 
displaced Cs ion has 3 nearest neighbours and a bond 
length of 3.8 A which from simple electrostatics gives an 
increase of 0.62 eV. The remaining energy changes are 
likely due to next nearest neighbours and subtle changes 
in the magnetism of the system. 

The two different stacking sequences proved to be ex¬ 
tremely close in energy. The total energy difference be¬ 
tween the two models was 1.5 meV. This is extremely 
small and only a little larger than the convergence of the 
calculation. This difference in energy equates to an or¬ 
dering temperature of 20A, implying that above this 
temperature there will be no correlations between sub¬ 
sequent layers, in agreement with the simulations of the 
Bragg rods above. At such a low ordering temperature, 
the large Cs ions are unlikely to be mobile and therefore 
the inter-layer disorder is likely to freeze in, explaining 
why no ordered version of this phase has been observed 
at low temperature. 


VII. DISCUSSION 

There have been a large number of studies on the 
phase-separation in this class of systems, with mi¬ 
croscopic phase separation being identified early on 
as a potential solution to the observed coexistence 
of superconductivity and magnetic ordering in the 
systenP. Multiple phases have been observed in all 
the variatio ns of A^F 62-^8 02 witlWA = A, Rb, Cs) 
using T EivPlp STlvp EZl^ NMEP and diffraction 

In each case, the system is 
found to form two, spatially separated phases, where the 
first, dominant, phase is identified as the a/5 x ^/b an¬ 
tiferromagnetic phase. The full structure of the second 
phase is, until now, unidentified. The size of these re¬ 
gions is smaller than that reported here, though STM 
and TEM techniques measure at a smaller length scale, 
implying that there is a patterning within the struc¬ 
tural domains that we have measured. A refinement of 
the second phase was performed from powder neutron 
diffraction measurements by Pomjakushin et al giving a 
composition in M/mmm of 61 ( 5 )^ 62 . 2 (i) 5 'e 2 , which 
is consistent with agreement with our refinements and 
modeP^. 

Synchrotron x-ray analysis by Ricci et al established 
the presence of a a/2 x ^/2 periodic superstructure^ and 
Bosak et al. observed the ID rods and a monoclinic 
distortiorP^ in this phase. These rods are consistent 
with our measurements, however we did not observe a 
similar monoclinic distortion on SXD due to the coarse 
wavevector resolution of the detectors. During our own 
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synchrotron measurement we did observe a splitting of 
the a/2 X a/2 peaks consistent with the monoclinic dis¬ 
tortion, however for our surface scans we integrated over 
the different components. It was also noted in the report 
by Bosak et al that Phase 2 is I centred, precluding dis¬ 
order in the FeSe layer stacking and the postulation was 
made for uncorrelated Cs layer 


In the case of RhxFe 2 -ySe 2 , polarised neutron scatter¬ 
ing has been used to determine that the \/2 x \/2 diffuse 
rods have no magnetic scatterin^^. NMR data by Texier 
et al. on RbxFe 2 -ySe 2 indicates the presence of a ma¬ 
jority magnetic phase in the presence of a minority, Fe 
vacancy free superconducting phas^^^. By probing the 
superconducting band gap in KxFe 2 -ySe 2 with a STM, 
Li et al determined that Fe vacancies are destructive to 
superconductivity in the systenP^, in agreement with the 
filled Fe layer in our model. Further analysis of STM im¬ 
ages by Cai et al. has suggested that the \/2 x \/2 phase 
contains Se charge orderin^^. On its own however, such 
charge ordering could not account for the Bragg rods ob¬ 
served here as such a signal would be vanishingly weak 
using neutron diffraction. It is possible, however, that 
charge ordering is driven by in-plane ordering of the al¬ 
kaline ions. The previous literature is therefore in good 
agreement with our Cs vacancy model for the secondary, 
superconducting phase. 


A recent soft chemistry study by Sun et al. has con¬ 
cluded that the highest superconducting temperatures in 
FeSe systems are achieved when the iron oxidation state 
is reduced below +2 and when the iron vacancy concen¬ 
tration is lovP^. This is consistent with the lack of va¬ 
cancies and reduced oxidation state implicit in our model 
with nominal concentration Cso,^Fe 2 Se 2 . 


VIII. CONCLUSION 

Single crystal neutron Lane diffraction and syn¬ 
chrotron x-ray diffraction have been combined to deter¬ 
mine the separation of phases in Cso,sFei,QSe 2 . Refine¬ 
ments of the two overlapping phases indicate that one 
phase comprises Fe vacancies and the other Cs vacancies. 
The two phases have different superlattices, with the Fe 
vacancy phase having a ^/E x \/5 superlattice and the 
Cs vacancy phase exhibiting a ^/2 x ^/2 superlattice with 
Bragg rods along 1. The two phases were determined to 
be spatially separated by x-ray diffraction measurements 
across the surface of the sample. The first phase can be 
attributed to the Fe vacancy phase well studied in the 
literature, however the Cs vacancy phase can now be de¬ 
termined as a 50% Cs phase with Cs ions sitting at every 
other allowed site. First principles calculations were able 
to determine that while the 2D layers are stable, there 
are no correlations between the inter-layer ordering of 
the Cs ions, giving rise to ID Bragg rods in reciprocal 
space. The observation of this second, spatially sepa¬ 
rated phase in the presence of a primary phase known 
to be highly magnetically ordered suggests that this new 
phase is contributing to the observed superconductivity 
of the system. 
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